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ABSTRACT 


A combination of several measurement techniques was used to investigate the 
dielectric properties of 80 rock samples in the microwave region. The real part of the 

dielectric constant, e\ was measured in 0.1 GHz steps from 0.5 GHz to 1 8 GHz, and the 

imaginary part, e", was measured at five frequencies extending between 1 .6 GHz and 
1 6 GHz. In addition to the dielectric measurements, the bulk density was measured for 
all the samples and the bulk chemical composition was determined for 56 of the 

samples. This study shows that e* is frequency-independent over the 0.5-1 8 GHz for 

all rock samples, and that the bulk density p accounts for about 50% of the observed 

variance of e'. For individual rock types (by genesis), about 90% of the observed 
variance may be explained by the combination of density and the fractional contents of 

SiC> 2 , Fe 203 , MgO, and TiC> 2 . For the loss factor e", it was not possible to establish 

statistically significant relationships between it and the measured properties of the rock 
samples (density and chemical composition). 
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1. INTRODUCTION 


The purpose of this study is to investigate the microwave dielectric properties of 
igneous and sedimentary rocks, in support of radar investigations of the Earth's 
geology and of the proposed Mars Orbiting Radar and Radiometer (MORAR) mission 
to Mars. 

Several studies have been reported in the literature on the dielectric properties of 
rocks [1-6], but in most of these studies the reported experimental measurements had 
been made either at MHz or lower frequencies, or at one or very few microwave 
frequencies. Thus, no continuous microwave spectra of the relative dielectric constant 

e have been reported to date. Furthermore, the majority of the reported data for the 

dielectric loss factor e" is of questionable accuracy. This is because e" of most rocks is 
between 0.01 and 0.1 , and most dielectric measurement techniques do not have the 
accuracy required for measuring values that small. 

The relative dielectric constant e of a material is defined as 
e * e’ - j e" , 

where the real part e' is the permittivity of the material (relative to that of free space) 

and the imaginary part e" is its dielectric loss factor (also relative to eq of free space). 
This study focuses on the spectral region extending from 0.5 GHz to 18 GHz. A 

combination of several measurement techniques was used to measure e over this 
frequency range. It included two probe techniques for measuring e’ in steps of 0.1 GHz 
from 0.5 GHz to 1 8 GHz, and a resonant cavity perturbation technique for measuring e" 
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at five frequencies extending from 1 .6 GHz to 16 GHz. As discussed later, because 
these cavity measurements are very time-consuming, it was not possible to make the 
measurements at more than five frequencies (within the constraints of available 
resources) without compromising measurement accuracy. 

The dielectric data reported in this study were generated from measurements 
performed for 80 rock samples. Each data point represents the average of several 
measurements corresponding to spatially different parts of the rock sample. The 
variability among measurements made for a given rock sample is an indicator of the 
sample's spatial inhomogeneity. Such variations may be due to density variations or 
variations in chemical composition among mineral constituents. In addition to 
measuring the dielectric behavior of each sample, its density and bulk chemical 
composition were measured and documented also. 

This report provides listings of the measured data, analyses of the associated 

measurement accuracies, and analyses relating e' and e" to the density p of the 
measured samples and to their chemical contents. 

2. DIELECTRIC PROBE MEASUREMENT TECHNIQUE 

The permittivity data reported in this study are based on measurements of the 
complex reflection coefficient of a coaxial probe terminated in the material under test. 
Two techniques were used. The first one is based on a third-order equivalent circuit 
that can be used for measuring the dielectric constant of any rock sample across the 
full frequency range of interest (0.5 - 18 GHz). The second one is a simpler first-ordei 
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equivalent circuit, but its validity range is limited to frequencies below 10 GHz if £' is 


larger than 8. For all rock samples investigated in this study, e' was found to be 
approximately independent of frequency over the 0.5-1 8 GHz range. Because it is 

simpler to use and calibrate, the first-order technique was initially used to measure e' 

of a given sample, and if e’ was found to exceed 8 over the 0.5-10 GHz range, the 
sample was remeasured using the more exact third-order technique. Brief 
descriptions of these two techniques are given next. 


2.1 Third-Order Equivalent Circuit 

The dielectric probe system (Fig. 1 ) consists of a swept RF source, a network 
analyzer (HP 851 0A), and associated couplers and data processing instrumentation. 
Fig. 2(a) shows a cross section of the probe tip and the dimensions of two of the 
probes examined in this study. The operation of open-ended coaxial lines to measure 
the dielectric constant of unknown materials is well-documented in the literature [7]-[9]. 

The input reflection coefficient at the prober tip p is given by 


Z L- Z 0 

\ + Z 0 


y o~ y l 
y o + y l 


( 1 ) 


where Y = 1/Z, Zq is the line impedance, and Zl is the load impedance, which is 
governed by the geometry of the probe tip and the dielectric constant of the material it 
is in contact with or immersed in (for liquid materials). In general, an open-ended 
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Block diagram of the dielectric-probe measurement system. 
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Fig. 2. The dielectric probe; (a) cross-sectional views, and (b) equivalent circuit. 
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coaxial line may be described by an equivalent circuit of the form shown in Fig. 2(b). 
When placed in contact with a homogeneous material whose thickness is sufficient to 
simulate a slab of infinite electrical thickness, an open coaxial line has an admittance 

Y|_(a), e) given by 

Y l (co, e) * Y. (co) + Y e (co, e) , (2) 

where Yj (co) ■ jcoCj is the "internal" admittance corresponding to the fringing 
capacitance Cj that accounts for the fringing field in the Teflon region between the 
inner and outer conductors of the line. The "external" admittance Y e , which is a 

function of both co and the complex dielectric constant e of the material under test, 
consists of a frequency-dependent capacitor C(co, e) in parallel with a radiation 
conductance G(co, e) 

Y (co, e) - jcoC (co, e) + G (co, 0 . (3) 

The capacitor C(oo, e) represents the fringing field concentration in the dielectric 

medium (e) surrounding the probe tip, and the conductance G(co, e) represents the 
radiation into the dielectric medium. 

When the medium surrounding the probe tip is free space (i.e., an open-ended 
line), these two equivalent-circuit elements vary according to 
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(4) 


2 

C (co, £q) = Cq + Bo) 

4 

G (co, £q) — Aco , (5) 

where Co, B, and A are constants for a given probe-tip geometry. If the radial 
dimensions of the coaxial line (namely, r-\ and r2) are small compared to the 

wavelength X, computations using the expressions given in Marcuvitz [1 0] yield values 
for A and B that are sufficiently small that the external admittance may be 

approximated as Ye(co, eq) = jcoCg. If the dielectric constant of the medium surrounding 

the probe tip is not the free space value eq, however, the above simplification may lead 
to unacceptably large errors. Hence, in the general case we have 


Y e (co, £q) = jco(C 0 + Bo) 2 ) + Aco 4 . (6) 

According to the theorem developed by Deschamps [11], the input admittance of 

an antenna immersed in a medium of complex dielectric constant e is related to the 
input admittance in free space through 



(7) 


The above expression is for materials characterized by p = jiq. If we regard the 
open-ended coaxial line as an antenna and henceforth abbreviate the relative 

dielectric constant ratio eJeq as simply e, we can write the following expression for the 
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total input admittance of the probe when placed in contact with a material of relative 


dielectric constant e: 

3 2 4 2 5 

Y l (cd, e) = jtoCj + jwC Q e + jBco e + Aco e (8) 

With the line admittance Yq known, measurements of the amplitude and phase of p by 
the network analyzer system (Fig. 1 ) lead to a measurement of Y|_. The next step is to 

determine e from Yl. This is accomplished by 1 ) calibrating the measurement probe in 
order to establish the values of the constants Cj, Co, B, and A, and 2) developing an 

iterative program for finding a value for e that minimizes the error between the 
measured value of Y l and the value calculated from the expression on the right-hand 
side of (8). 

A Usable Frequency Range 

The radii r-j and r 2 of the coaxial line govern three important characteristics of the 
dielectric measurement system: 

1 ) The ratio r*|/r 2 determines the characteristic impedance Zq of the line. For 
50-ft Teflon-filled lines, this ratio is approximately 0.3. 

2) The difference 02 - r\ ) determines the cutoff wavelength of the TM modes 

[1 0, p. 74]; the cutoff wavelength of the TMqi modes is Xq =2{r2 - r-| ). Table 1 
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Table 1. Dimensions and cutoff wavelength X c for the TMqi mode for 
four standard-size coaxial cables. 


Cable 

Type 

^(mm) 

0.09" 

Teflon 

0.84 

0.14" 

Teflon 

1.50 

0.25" 

Teflon 

2.66 

0.35" 

Teflon 

3.62 


ri(mm) 

r 2 /r 1 

X c (mm) 

0.26 

3.28 

1.18 

0.46 

3.30 

2.13 

0.82 

3.22 

3.76 

1.12 

3.22 

5.07 
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provides a list of the dimensions and cutoff wavelengths of four standard coaxial lines 
used in this study. For a medium with a complex dielectric constant e ■ e' - je" , the 


wavelength in the medium Xg is related to the wavelength in free space Xq by 

,-1/2 


V x o 


§• (1 + sec 5) 


(9) 


where 5 = tan*1 (e"/e'). To avoid the propagation of TM modes, the condition X e < X c 

should be satisfied. Because most rocks are low-loss materials and s' seldom 

exceeds 10, all of the probes listed in Table 1 are appropriate for measuring e' of rocks 
at frequencies below 20 GHz. In this investigation 0.1 4-in probes were used 
exclusively. 


B. Calibration 

Calibration entails finding the values of the constants Cj, Co, B, and A of (8) for 
each probe used in this study. Under ideal circumstances, one needs to determine 
these constants only once and at only one frequency. The equivalent-circuit model, 
however, is only approximate; hence, it is necessary to determine these constants at 
each frequency that the probe is intended to be used. For example, it was found that 

the constant A varies approximately at 1/co, which means that the conductance term 

G(co) varies as co 3 , not co 4 . 

Each dielectric probe was calibrated by measuring the complex reflection 
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coefficient under four termination conditions: 1 ) short circuit, 2) open circuit, 3) probe 
immersed in distilled water, and 4) probe immersed in methanol. Distilled water and 
methanol were used because their dispersion spectra are well known [12], [13]. 

2.2 First-Order Equivalent Circuit 

If the diameter of the coaxial probe is much smaller than the wavelength in the 
material under test, the equivalent-circuit admittance Y e simplifies to only one term, 

jooCoe, because the other two terms become negligibly small. For the 0.14-in probe 
used in this study, the condition 

f(GHz) <; 50 /V? . ( 10 ) 

must be satisfied in order for the first-order model to yield accurate results. This 
condition was found by comparing measurements made with this technique to 
measurements made using the more-exact technique described in the previous 
section. 

2.2.1 Reflection Measurement Technique 

For the first-order equivalent circuit, the admittance, 

Y l (co, e) = j co Cj + j to C Q e , (11) 

can be determined by measuring the reflection coefficient p. The constants Cj and 
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Cq can be determined by measuring Yl(co, e) for two materials with known e. With the 
constants known, e of an unknown material may be computed directly from 






( 12 ) 


by measuring p. The coaxial line is a standard 50-ohm line (i.e., Yq = 1/50). 


2.2.2 Group-Pelav Measurement Technique 

As an alternative to measuring p in order to determine e, a group-delay 
technique was developed which requires calibration against only one calibration 

material rather than two. For low loss materials with e" « e\ 


Y l = jatC. + e’ C Q ) 


and the reflection coefficient 


P = Ipl 0 


j<t> 


Y - Y 

_0 _L 

Y + Y 
0 L 


has a phase angle given by 


-1 


<t» = 2 cot [50 a) (C. + e' C n )] 

I u 


(13) 


(14) 


(15) 
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For the 0.1 41 -in probe, the constants Cj and Co are on the order of 0.02 picofarads. 
Consequently, the entire quantity inside the square brackets is much smaller than 1 if 

co£2tcx20GHz and e'^10. Hence, cot'^ ) may be expanded in a Taylor series 

cot 1 (x) = 7C / 2 - x + x^/3 - . . . , (16) 

and if we retain only the first two terms, we have 

4> s jc - 1 00 co (Cj + e'C 0 ) . (17) 

The group delay t is defined as the change in the phase of p as a function of 
frequency, 

x = = -100 (C. + e' CJ . (18) 

dco 1 0 

If the group delay is measured with the probe in air (with e' = 1 ) and not in contact with 
any other material, we get the reference group delay xg, 


x 0 = -100 (Cj + C Q ) 


The differential group delay is defined as 


Ax = x - x. 


= - 1 00 C Q (e' - 1) , 


(19) 


( 20 ) 
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from which we obtain the expression 


e' 


Ax 

100 C 0 ' 


(2i ; 


The constant Cq may be determined by measuring Ax for one material of known e'. 

The group delay t and xq can be measured directly by the HP 851 OA network 
analyzer. 

Comparison of results using the group-delay technique with results obtained 
using the more-exact reflection coefficient technique has led to the conclusion that the 
condition 


f (GHz) <; 30 /J7 , (22 

should be satisfied in order for the approximation made in going from (16) to (17) to be 
valid. 


2.3 Sample Preparation For Permittivity Measurements 

When using the coaxial probe to measure the permittivity of a solid material, the 
following two conditions must be satisfied (in order for the measurements to produce 
accurate results): 

(1 ) The thickness of the sample must be at least equal to the probe diameter. 
For the 0.1 4-in probe, this condition is satisfied if the thickness is greater than 4 mm. 
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(2) The surface of the sample in contact with the probe must be very smooth in 
order to insure good electrical contact. This was achieved by having each rock 
sample cut with a rock saw to obtain a flat surface and then the surface was smoothed 
using a table-top rotary sander. 

To avoid dielectric effects that may be caused by the possible presence of 
surficial water molecules on the sample, each sample was dried in an oven for 15 
minutes at 100°C prior to performing the dielectric measurements. It was found, 
however, that there was very little difference, if any, between the results obtained after 
drying the samples and those obtained on the basis of the measurements made prior 
to drying the samples. An entirely different conclusion was reached for the 

measurements of the dielectric loss factor e"; for some rocks, the values measured 
prior to drying the sample were as much as twice the values measured for the samples 
dry (Section 3.3). 

2.4 Measurement Accuracy and Precision 

The measurement accuracy of the probe technique was evaluated by comparing 
the permittivity measured by the probe with the permittivity of standard materials. The 
reference materials are homogeneous, thick blocks of solid materials, such as teflon, 
whose dielectric constants had been carefully measured using waveguide techniques. 

A typical comparison is shown in Fig. 3 for a material with e’ * 8.0. Based on this and 
other comparisons, we estimate the probe measurement accuracy to be better than 
±0.03 of the measured value. 
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Fig. 3. Comparison of the measured permittivity with the permittivity of a reference 
material. 




By accuracy, we refer to the absolute level of e', whereas by precision, we refer to 
the variability associated with the spatial inhomogeneity of the sample. For all rock 

samples, e' was measured by applying the probe at at-least 16 spatially different 
locations on the polished surface of the rock sample. In each case, we computed the 

mean value of e' , the associated standard deviation S, and the ratio S/ e’. For 79 of 

the 80 samples, the ratio S/ e' was found to be smaller than 0.12, and for most the ratio 
was smaller than 0.05. The one exception was rock sample PW-30 (rock samples are 

identified by an identification number as shown in Table 3) for which e' varied from 6.6 
to 13.2. The reason for the variability was visible on the rock's surface; Fig. 4(a) shows 
a sketch of the surface of sample PW-30, which is a two-tone material comprised of a 
light-tone background and darker-tone inclusions. The sketch in Fig. 4(b) shows the 
permittivities measured at six locations on the surface of the sample; for locations not 
including the inclusion material, the measured value was 6.6, whereas for those 
partially or totally covering the inclusion material, the values were higher. Because of 

this large spatial variability in e', the sample was cut and only the "homogeneous" 
background portion was used in the analysis. 

Figure 5 shows typical permittivity spectra for four rock samples, which exhibit no 
discernible dependence on frequency. This is characteristic of all samples measured 
in this study and is in agreement with previous conclusions reached by Olhoeft 

et al. [2]. Hence, in all forthcoming discussions and analyses, e' will be treated as 
frequency independent and will be represented by the average value measured over 
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(a) Sketch of rock surface (b) Measured permittivity 


Fig. 4. This figure illustrates the spatial variability of the permittivity of sample 

G-PW-30. The dark areas in (a) consist of high iron content inclusions, and 
the circles in (b) show the locations at which the probe measurements were 
made. 
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Permittivity 
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Fig. 5. Typical examples of measured permittivity spectra. 



the 0.5-18 GHz range. 


3. RESONANT-CAVITY PERTURBATION TECHNIQUE 

A resonant cavity is a closed volume. Metal couplers protruding slightly into the 
cavity volume are used to measure the resonance characteristics of the cavity. The 
diagrams in Fig. 6 show a cylindrical cavity with two magnetic loop couplers protruding 
slightly into the cavity volume on the inside walls and connected to SMA connectors 
on the outside walls at a height midway between the top side (the lid) and the cavity 
floor. Figure 7 shows the measurement system. 

With the cavity empty, if one were to connect a signal generator (HP 8350B in Fig. 
7) to one of the connectors and a network analyzer to the other and then sweep the 
generator frequency across the resonance region of the cavity, the output power would 
be a Gaussian-like function of frequency (Fig. 8). This power spectrum is 

characterized by fg, the frequency at which the power is a maximum, and by Qg, the 
quality factor, 



(23) 


where Af is the half-power width of the power spectrum. If we insert a dielectric 
material into the cavity, the spectrum will change in two ways: 1 ) the resonant 
frequency decreases to a lower value, which we shall call f s , and 2) the quality factor 

decreases to a lower value Q s . 
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Fig. 6. Diagrams showing the coupling arrangements used with the cylindrical cavity. 





HP8510 SYSTEM BUS 



Fig. 7. Block diagram of the system used to measure the cavity transmission spectri 
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Sample 



Resonant Cavity 
(a) Cavity with Sample 



(b) Transmission Spectra 

Fig. 8. When a sample is inserted in the cavity, its transmission spectrum changes: 
the resonant frequency shifts from fg to f s and the spectrum becomes 
broader. 
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In order to maintain Q s large (i.e., maintain a resonant-like spectrum), the 

volume of the material inserted into the cavity must be kept small relative to the cavity 
volume. When this is the case, the resonant-cavity perturbation technique [14] may be 

used to determine e' and e" of a dielectric material from measurement of fo, f s . Qo. 
and Q s . 

For a cylindrical cavity with radius a and height d<2a oscillating in the TMqio 
mode and containing a needle-shaped dielectric material oriented along the vertical 
axis of the cylinder, the shift in the resonant frequency is 


S ° - - 1.855 V(e* - 1) (24) 

T s 

if the volume fraction V is small. If the material has a dielectric loss factor e", it can be 
shown [1 4, p. 373] that 


e" 


w z H 
2 





Solving (24) for e' we get 


(25) 


e' = 1 + 


f 0- f s 
1.855 f V ‘ 


( 26 ) 


These expressions are valid only if V is very small, and (24) is valid only if the 
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dielectric material is approximately needle shaped and oriented vertically. One of the 


major problems associated with using this method to determine e' and e" is the need to 

know V very accurately (e" depends on (e* - 1) which, in turn, depends directly on 1/V). 

In our case, however, we did not need to know V exactly because we already know e’ 

from the probe measurements discussed in the previous section. Hence, e M could be 
determined from (25) without the need to measure V. This procedure of using 

dielectric probes to measure s' and resonant cavities to measure e" proved extremely 
effective because the errors associated with the handling and the measuring of the 
weight and volume of very small rock samples were intolerably high. As will be 
discussed below, a desirable value for V is about 0.5 percent. For a cavity volume of 
2.5 cm^ (which was the volume of one of the cavities used in this study), V would have 
to be about 1.25 x 10' 2 cm 3 and the corresponding weight would be about 31 mg (for 
a typical density of 2.5 g/cm 3 ). 


3.1 Measurement Accuracy a nd Precision 

By way of evaluating the measurement technique as well as establishing the 
range of validity of (25) as a function of V, we conducted a carefully designed 

experiment in which e" of plexiglass was measured as a function of V for values of V 
extending from 10* 3 percent to 10.1 percent. We chose plexiglass because its 

complex dielectric constant is well known (e - 2.55 * j 0.0165) and its dielectric loss 
factor is small. The results of the experiment are shown in Fig. 9. The measurement 
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technique predicated the correct value for e" within a rms error of 0.001 for the range 
0.01 percent £ V £ 1 percent, and with a slightly larger error up to 5 percent. A detailed 
analysis of the errors associated with the measurements of the quantities fo, f s , Qq, 

and Q s led to the conclusion that the optimum range of V is between 0.5 percent and 

1 percent, and that if V is in this range the minimum measurable value of e N is around 
0 . 002 . 

The data tabulated in Table 3 are each an average of measurements conducted 
for five small samples of the parent sample. 

3.2 Cavltv Characteristics 

Five cylindrical cavities were used in this study, with center frequencies ranging 
from 1 .6 GHz to 1 6 GHz. Table 2 provides a listing of their pertinent characteristics. Of 
particular note is the cavity volume, ranging from approximately 1000 cm 3 for the 
1 .6 GHz cavity to only 1 cm 3 for the 1 6 GHz cavity. 


3.3 Effects of Surflctal Water 

The data provided in Section 4 are based on measurements conducted after 
drying each rock sample in an oven at 100°C for 24 hours. The samples were dried to 

avoid the effects of surficial water on e". For most of the samples, no discernible 

difference in e" was observed between measurements made before and after drying 
the sample (Fig. 10(a)), while for some the difference was quite significant (Fig. 10(b)). 
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Table 2. Characteristics of resonant cavities 


Cavity Frequency (GHz) 

1.6 

5 

7.8 

11.4 

16 

Resonant Frequency (GHz) 

1.64 

4.98 

7.87 

11.44 

15.93 

Qq (empty cavity) 

3000 

1150 

480 

270 

55 

Height (cm) 

6.34 

2.06 

1.25 

0.82 

0.65 

Diameter (cm) 

13.93 

4.57 

2.91 

1.97 

1.43 

Volume of cavity (cm3) 

966.5 

33.75 

8.3 

2.5 

1.05 
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Measured densities and dielectric constants of rocks. 
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For a few of the measured samples, e" was observed to exhibit no discernible 


dependence on frequency. For most samples, however, e" decreased with increasing 
frequency over the 1 .6 GHz - 1 6 GHz range. Typical examples of these two types of 
spectra are shown in Fig. 1 1 . 

4. MEASURED DATA 

The measured dielectric data is given in Table 3. The entries include 1) rock 
type, 2) rock #, which designates the source (G = NASA/GSFC, J = NASA/JPL, and 

E a ERIM) and associated numbers, 3) density, 4) e' (average value over the 

0.5-18 GHz range), 5) S/e', the standard derivation-to-mean ratio of the measured 

value of e\ and 6-10) are entries fore" at 1.6, 5.0, 7.8, 1 1.4, and 16.0 GHz. 

Table 4 lists the bulk chemical properties of the rock samples from x-ray 
fluorescence. The properties of samples with a G designation were measured by or 
for NASA/GSFC, and the properties of those with a J designation were measured by 
the Geology Department at the University of Michigan. To date 56 of the 80 rock 
samples have been analyzed for mineral composition. 

5. ANALYSIS OF PERMITTIVITY DATA 
5.1 Distribution of Measured Data 

Among the 80 rock samples, the measured value of e' ranged between 2.5 and 
8.3. These values are presented in horizontal bar-chart format in Fig. 1 2, and a similar 
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Dielectric Loss Factor 


0.3 


-0- G-2 

J-7-21-84L6 



Frequency (GHz) 


Fig. 1 1 . Typical examples of the measured spectra of the dielectric loss factor e". 
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Table 4. Mineral composition of the rock samples. 
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presentation is given in Fig. 13 for the rock groups identified in Tables 3 and 4. 


5.2 Dependence on Density 

According to previous studies [1-6], the density p(g/cm3) j s the single most 
important parameter governing the magnitude of e\ One of the commonly used 

formulas relating e* to p is based on a simple model relating of a powder 

material of density pp to that of the parent (solid) rock material of permittivity e and 
density p through 


% 1/p 1/p a 

V - « p 


= (e) 


1/p 


(27) 


where e a = 1 and p a = 1 are the permittivity and density of air. Campbell and Ulrichs 

[1] conducted measurements for a large number of powdered rocks, all at a density 
3 

pp= 1 g/cm , and found that varied over the narrow range between 1 .9 and 2.1 for 

most of the 25 different types of powdered rocks measured and that the mean value is 

3 

around 2.0. Upon setting p p * Ig/cm and e’ p * 2 in (27), we get 


= 2 p . 


(28) 


This result is in close agreement with the formula used by Olhoeft and Strangway [4], 


e* = (1.93 ± 0.1 7) p , 


(29) 


in their analysis of moon rocks. For the data measured in the present study, an 
equation of the form e' = A p was used to fit the data and the value A = 1 .96 was found 
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o 
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Fig. 13. Range of e' for individual rock types. 



( 30 ) 


to give the minimum mean-square-error. The function 
e' = 1.96 p , 

is shown in Fig. 14, together with the measured data. The linear correlation coefficient 

between the values predicted by (30) and the measured values of e' is R = 0.72. The 

data scatter about the regression curve is attributed to the dependence of e‘ on the 
mineral composition of the rocks. 

Also, the data was used to generate first-, second-, and third-order polynomial 
regressions relating (e* - 1 ) to p. The results are: 

e’-1 - 1.86 p , (31a) 

e'-1 = 0.61 p + 0.49 p 2 , (31b) 

e‘-1 = 1.0 p- 0.11 p 3 , (31c) 

and in all cases the linear correlation coefficient between the measured value of e' 
and the value predicted by any of the above equations was R = 0.73. Thus, 
statistically speaking, the model given by (30) offers a fit to the data comparable to that 
provided by the simple linear model given by (31a), which is shown in Fig. 15. 

The variation of e’ with p is shown in Fig. 1 6 for individual rock types. The 
carbonates exhibit the narrowest density range, followed by the igneous plutonic 
silicates, and then by the igneous volcanic silicates and the sedimentary silicates. The 

igneous volcanic silicates have the strongest slope for s' versus p. 
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p 

Measured and Predicted (1.96 ) vs. Density 




Fig. 14. Measured permittivity and predicted permittivity = 1.96 p , both 
plotted against density in (a) and against each other in (b). 
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Predicted Permittivity 


Measured and Predicted (1 + 1.86p) vs. Density 



Density (g/cm3) 

(a) e’ versus density 


Predicted (1 + 1.86p) vs. Measured Permittivity 

9 
8 
7 
6 
5 
4 
3 
2 
1 
0 

01 23456789 

Measured Permittivity 

(b) e' versus e’ 
p m 

Fig. 15. Measured permittivity e' m and predicted permittivity e'p= 1 + 1.86p, both 
plotted against density in (a) and against each other in (b). 
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Permittivity versus Density 



Density (g/cm3) 


Fig. 1 6. Permittivity versus density for individual rock types. 
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5.3 Dependence on Mineral Composition 


Stepwise multivariate regression analysis was performed (using BMDP) to 

evaluate the statistical correlation between the measured e’ of the rock samples and 
their bulk chemical composition. The analysis was performed for individual types of 
rocks as well as for combinations of rock types. 

5.3.1 Igneous Volcanic Silicates 

For the 20 silicate samples classified as igneous volcanic rocks, the stepwise 
regression analysis selected density as the most important variable. According to 
Tabfe 5, density accounts for 76% of the total variance (see the R-squared entry in the 
table), the combination of density and Si02 content accounts for 88% of the variance, 
and the other elements are each responsible for small incremental improvements in 
the variance. With density and Si02 content, we obtain a linear correlation coefficient 

of 0.94 between the measured value of e' and the value computed using the linear 
regression equation given in Fig. 17(a). 

5.3.2 Igneous Plutonic Silicates 

A similar analysis conducted for the igneous plutonic silicate samples gives the 
results tabulated in Table 6. Due to the high correlation between iron content and 
density, density was not selected directly as a significant variate. Figure 1 7(b) shows 

the results of regressing e' against a linear equation containing total iron oxide 
content (Fe 203 T), Na20, and Si02, as parameters. The plutonic results in Table 6 
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Table 5. Stepwise multivariate regression analysis for e* for igneous 
volcanic silicate. (N = 20). 
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(a) Igneous Volcanics 


Fig. 1 7. Computed permittivity versus measured permittivity for (a) igneous volcanics 
and (b) igeneous plutonics. 
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Table 6. Stepwise multivariate regression analysis lor e’ for igneous 
Plutonic silicates. (N = 14). 
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£p = 2.5657 > 0.2443(Fe203T) + 0.1753(Na2O) + 0.0255(SiO2) 
(b) Igneous Plutonics 
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show that the mass percentages of a greater number of cations are significant in 

determination of e' than was found to be the case for the volcanics in Table 5. This 
implies that mineralogy plays an important, through possibly secondary role to density, 

in determining e'. 

5.3.3 Igneous Silicates 

If we combine all 34 igneous silicates into a single class, application of stepwise 
multivariate regression analysis provides the results given in Table 7 and Fig. 18. 

5.3.4 Sedimentary Silicates 

For the 20 sedimentary silicate samples, MgO content was selected by the stepwise 

regression program as the most important variable, followed by Si02. AI 2 . 03 , Ti 02 , and 
density. Results of the regression analysis are given in Table 8 and Fig. 19. Table 8a 
and Fig. 1 9a shows the results obtained when density and oxides are included in the 
analysis. These results are changed dramatically by inclusion of LOI (loss on ignition to 
1 ,000°C) in the analysis as shown by Table 8b and Fig. 1 9b. Loss on ignition is related 
to the vaporization of volatiles, including chemically bound water, during sample 
preparation for x-ray fluorescence studies. It is probable that the large amount of 

variance in e’ explained by LOI (56%) for the sedimentary silicates is related to the 
presence of chemically bound water in the sediments. It is interesting to note that 

inclusion of LOI in the analyses of e’ for the igneous silicates did not statistically alter the 
results. 
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e p - 2.0387 + 0.1 154(Fe203T) + 2.0947(Density) - 0.0591 (Mg O) 


Fig. 1 8. Computed permittivity versus measured permittivity for all igneous silicates 
combined. 
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(b) includes LOI in analysis 







































































e'p = 6.2882 + 1.6235(MgO) -0.046 1 (Si0 2 ) + 0.3820(AI 2 0 3 ) 
- 2.281 9(Ti0 2 ) 


(a) Without including LOI (loss on ignition) in analysis 


Fig. 1 9. Computed permittivity versus measured permittivity for sedimentary silicates 
(a) without including LOI (loss on ignition) in analysis and (b) including LOI 
in analysis. 
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£p » -1 .7925 + 0.31 73(LOI) - 0.1 224(CaO) + 0.4442(Na 2 O) 
+ 2.0452(density) 


(b) Including LOI in analysis 
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5.3.5 All Silicates 


Of the total of 54 silicate samples, 34 were igneous and 20 were sedimentary. 

The combined analysis in Table 9a shows density, aluminum and iron content to be 
the most significant variates. Not suprisingly, a significant portion of the variance 
within each of the silicate subgroups is explained by one or more of these variates. 

Figure 20 shows a plot of a linear regression equation relating e’ to density and bulk 

chemical properties for the silicate samples. It is apparent that e' is dominated by the 
density of the rock and is secondarily influenced by the chemistry of the sample as a 
consequence of its mineralogy. Inclusion of LOI in this analysis does not significantly 
alter these results as shown by Table 9b. 

6. ANALYSIS OF DIELECTRIC-LOSS DATA 
6.1 Frequence Variation 

The dielectric loss factor e" was measured for 72 rock samples at five frequencies 

extending between 1 .6 GHz and 16 GHz. Figure 21 (a) shows plots of e" versus 

frequency for four rock groups. Each data point represents the average value of e" at 
that frequency for all rock samples belonging to that group. An overall-average plot for 

all 72 samples is shown in Fig. 21 (b). The plots in Fig. 21 indicate that e" decreases 
with increasing frequency between 1.6 GHz and 5 GHz and then it levels off at higher 
frequencies. Among the rock groups shown, the carborates exhibit the lowest loss and 
the igneous volcanics exhibit the highest loss. 
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Table 9. Stepwise multivariate regression analysis lor e’ lor all 

silicate rocks (a) analysis excludes LOI f and (b) analysis 
includes LOI. (N = 54). 
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W analysis includes LOI. 
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e' » -1 .9468 + 2.4689(Density) + 0.0786(AI203) + 0.0438(Fe203T) + 
P 0.0441 (CaO) - 0.3575(TiO2) 


Fig. 20. Predicted versus measured permittivity for all silicate rocks. 
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Predicted Dielectric Loss Factor vs. Frequency (GHz) 
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Frequency (GHz) 

(a) Individual Rock Types 

Fig. 21 . Spectral variation of e H for (a) individual rock types and (b) the average for all 
rock samples. 
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0.04 H 



The loss factor e" may be modeled as the sum of a conductive component e" and 

c 

a frequency-independent residual component e", 
e"(f) = e" + e“ (f) 


= P 1 + Pg/f ( 32 ) 

where o is the conductivity, and P 1 and P g are abbreviations for ej and c/2ne Q , 

respectively. For each of the measured samples, the values of Pi and P 2 were 
determined by fitting the measured data to a linear function of the form given by (32). 
The frequency was expressed in GHz, which makes P 2 have units of GHz* 1 . 

Figure 22 presents the data measured at 1 .6 GHz and 1 6 GHz in bar-chart format, 

arranged according to the magnitude of e" at 1 .6 GHz, starting with the largest at the 
bottom of the left-hand chart and ending with the smallest value at the top of the 
right-hand chart. Note that the scale is different for the two charts. At 1 .6 GHz, the 

magnitude of e" extends between a high of 0.24 and a low of less than 0.002. 

The constant P 2 of a given sample is proportional to its conductivity a. The 
values of Pi and P 2 determined by the 72 rock samples are arranged in bar-chart 
format in Fig. 23 according to the magnitude of P 2 . 
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Fig. 23. Bar chart presentation of the constants Pi and P 2 


6.2 Dependence on Density 


Unlike e' which exhibits a strong dependence on bulk density, the loss factor e" 
does not appear to be correlated with p(Fig. 24). The same conclusion applies to Pi 
and P2 (Fig. 25). 

6.3 Dependence on Mineral Composition 

Because e" is a function of frequency, it was decided to apply the stepwise linear 

regression analysis on Pi and P2 individually, rather than on e". The results are 
tabulated in Tables 10-14 and in Figs. 26-30. Among the 1 1 variables examined, 

Fez O3 was found to be the most important single variate for the igneous rocks (Table 
1 2), but that was not the case for the sedimentary rocks (Table 1 3). 

7. CONCLUSIONS 

The analysis of microwave dielectric measurements conducted for 80 rock 
samples has led to the following conclusions: 

(1 ) The permittivity e' is frequency-independent over the 0.5-1 8 GHz range. 

(2) The dielectric loss factor e" exhibits a frequency dependence of the form 

e" = Pi + P2 / f. For most samples, the second term is significant for f < 5 GHz, and 
may be ignored above 5 GHz. 
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Fig. 25. The constants Pi and ?2 versus bulk density. 






Predicted Dielectric Lass Factor Predicted Dielectric Loss Factor 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 


Measured Dielectric Loss Factor 
(a) 1.6 GHz 



0.00 0.05 0.10 0.15 0.20 

Measured Dielectric Loss Factor 
(b) 16 GHz 

e F> * P 1 +P 2 /f 

P 1 * 0.0341 + 0.0302(Fe203T) - 0.01 18(MgO) - 0.0090(CaO) 
- 0.2523(MnO) 

P 2 = - 0.0689 + 0.003(A12o©)- 0.0 117(Na20) + 0.0054(K 2 0) 
+ 0.0297(Density) 


Fig. 26. Predicted versus measured dielectric loss factor of igneous volcanic rocks. 
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Predicted Dielectric Loss Factor Predicted Dielectric Loss Factor 



Measured Predicted Dielectric Loss Factor 
(a) 1 .6 GHz 



Measured Dielectric Loss Factor 
(b) 16 GHz 

e" - P 1 + P 2 /f 

P 1 - 0.4559 + 0.223(MnO) - 0.1 724(Density) + 0.01 1 4(Fe203T) 

- 0.001 5(AI203) 

P 2 - 0.4685 + 0.0208(Fe203T) - 0.1789(Density) - 0.1734(P2O5) 

- 0.0024(MgO) 


Fig. 27. Predicted versus measured dielectric loss factor of igneous plutonic rocks. 
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Measured Dielectric Loss Factor 


(a) 1.6 GHz 



Measured Dielectric Loss Factor 
(b) 16 GHz 

e" * p i + P 2 /f 

P 1 = 0. 1 38 + 0.01 87(AI203) - 0.021 6(K20) -0.11 69(Ti02) 

- 0.001 4(Si02) - 0.0230(Density) 

P 2 - - 0.062 + 0.0392(Na20) - 0.01 66(K20) - 0.2706(P205) 

+ 0.0028(Fe203T) 

Fig. 28. Predicted versus measured dielectric loss factor for all igneous rocks. 
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Predicted Dielectric Loss Factor Predicted Dielectric Loss Factor 



Measured Dielectric Loss Factor 
(a) 1.6 GHz 



(b) 16 GHz 


e" * p i + P 2 /f 

P 1 = 0.1385 + 0.0131 (Fe203) - 0.0449(Density) - 0.0023(MgO) 

- 0.0301 (Ti02) 

P 2 * 0.0313 + 0.0154(Fe203) - 0.0082(CaO) - 0.0037(MgO) 

- 0.071 2(P205) 


Fig. 29. Predicted versus measured dielectric loss factor for sedimentary silicates. 
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(a) 1 .6 GHz 



(b) 16 GHz 

e" = Pi + P2/f 

Pi = 0.2224 -0.001 9(Si02)- 0.0040(K20) - 0.0026(MgO) - 
0.0062(Na20) - 0.0180(MnO) 

P 2 = -0.0381 + 0.0027(Fe203T) - 0.0772(P205) + 0.0037(AI203) - 
(0.0053) + 0.0007(Ai02) 

Fig. 30. Predicted versus measured dielectric loss factor for all silicate rocks. 
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Table 10. Stepwise multivariate regression analysis for e" for Igneous 
volcanic silicates. (N = 20). 
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Table 1 1 . Stepwise multivariate regression analysis for e M for igneous 

Plutonic silicates. (N = 14). 
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Table 1 2. Stepwise multivariate regression analysis for e" for all 
igneous silicate samples. (N = 20). 
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Table 14. Stepwide multivariate regression analysis for e" for all 
silicate samples. ( N = 54). 























































































































(3) The bulk density p accounts for about 50% of the variance of e*. For 

individual rock types (by genesis), about 90% of the observed variance of e’ may be 
explained by the combination of density, Sj02 content, Fe203 content, M g 0 content, 
and Tj O2 content. 

(4) The loss factor e" appears to be statistically uncorrelated with density. 

(5) Although multivariate regression equations were generated to relate e" to 
bulk sample chemistry, the results are not considered very significant (statistically) 

because the correlation coefficient between the measured and computed values of e n 
is poor. 

(6) Additional tests are needed to determine the dependence of e" on mineral 
composition. 
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